On the emergence of very long time fluctuations and 1// noise in ideal flows 
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This study shows the connection between three previously observed but seemingly unrelated 
phenomena in hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulent flows, involving 
the emergence of fluctuations occurring on very long time scales: the low-frequency 1/f noise in 
the power frequency spectrum, the delayed ergodicity of complex valued amplitude fluctuations in 
wavenumber space, and the spontaneous flippings or reversals of large scale fields. Direct numerical 
simulations of ideal MHD and HD are employed in three space dimensions, at low resolution, for 
long periods of time, and with high accuracy to study several cases: Different geometries, presence 
of rotation and/or a uniform magnetic field, and different values of the associated conserved global 
quantities. It is conjectured that the origin of all these long-time phenomena is rooted in the 
interaction of the longest wavelength fluctuations available to the system with fluctuations at much 
smaller scales. The strength of this non-local interaction is controlled either by the existence of 
conserved global quantities with a back-transfer in Fourier space, or by the presence of a slow 
manifold in the dynamics. 

PACS numbers: 47.27.-i,47.27.Sd,52.30.Cv, 52.35.Ra 



I. INTRODUCTION 



A characteristic of turbulence is the dynamical involve- 
ment of fluctuations over a broad range of length scales. 
At the largest scales the general expectation is that of 
non-universal behavior and influence by boundary condi- 
tions and/or driving. At intermediate scales, statistical 
behavior can be obtained, and at the smallest scales, be- 
havior is determined by some dissipation mechanism. As- 
sociated with these length scales, there are corresponding 
ranges of time scales. The characteristic nonlinear time 
scale, assuming local interactions in scale, is obtained di- 
mcnsionally as the ratio l/ui, where I is a length and 
ui is the corresponding typical velocity at /. Suppose we 
choose a particular length L and its associated time scale 
L/ul as the unit of time measurement. Frequently L will 
be the largest length scale in the system, or perhaps the 
size of the energy-containing eddies. In these units the 
range of locally computed time scales / /ui will typically 
extend from unity (or times of order one) to smaller time 
scales at the smaller length scales. However, there are 
some situations in which there is an emergence of much 
longer time scales, which cannot be associated with this 
kind of an estimate that is local in scale. Some of these 
issues have previously been identified in numerical simu- 
lations [l| , for example in the case of a driven-dissipative 
three-dimensional (3D) magnetohydrodynamic (MHD) 
system in presence of a strong background magnetic field, 
and in two-dimensional (2D) hydrodynamic (HD) and 
MHD driven-dissipative systems. In those three cases, 
the presence of long time fluctuations is connected with 
the appearance of a 1// noise type frequency power spec- 
trum at very low frequencies << u^/L. In many situa- 
tions SL 1/f noise spectrum arises from the existence of 



multiple correlation times in the system (see 0)- Such 
1/f spectra associated with long time scale fluctuations 
have been observed in a wide range of situations 0,3 but 
the presence of 1// noise in turbulent flows was reported 
only recently Dynamo simulations of the generation 
of magnetic fields also exhibit 1/f spectra Q. Low fre- 
quency dynamics have also been reported in laboratory 
experiments of HD and MHD fiows Closely re- 

lated are observations in the solar wind of 1// behavior 
in the spacecraft frame magnetic and density spectra at 
low frequencies @, Q . However in the solar wind case the 
1/f signal appears to trace back to the corona or even 
to the photosphere [l^l- It is therefore a consequence 
of long time scale variability of the solar wind sources, 
rather than a property that emerges due to solar wind 
dynamics itself. 

It is conjectured in [l| that the long time fluctuations 
in the turbulent MHD flows arise from the non-local cou- 
plings between the longest length scale in the system (the 
k = 1 mode in Fourier space) with the smaller scales, 
from the inertial to the dissipation range (fc S> 1). Such 
non-local interactions are known to be stronger for MHD 
than for HD flows ■ It remains however an open is- 

sue as to whether the existence of dissipation, the effects 
of the driving, or the particular choice of boundaries, are 
influential in determining the character or presence of the 
1/f noise that is generated. The present study attempts 
to show that the presence of a 1// noise spectral regime 
is specifically with the structure of the equations and the 
non-linear couplings. We accomplish this by employing 
a series of numerical simulations of ideal flows in a vari- 
ety of situations. We emphasize that these simulations 
are neither dissipative, nor driven. Furthermore, the 1/f 
regime will be characterized by showing its existence (or 
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absence), and its range of frequencies, for several different 
kinds of ideal flows and geometries. We will argue that 
a key feature in each system is the existence of either 
ideal global invariants, or quasi-invariants (to be defined 
below), which provide constraints on behavior of a small 
number of important degrees of freedom in the system. 

We will further argue below that a related issue con- 
cerning ideal flows is the phenomenon of delayed ergod- 
icity , also identified as broken ergodicity [ij, [T^ . In 
this phenomenon, the ideal flow is described as a statis- 
tical mechanical system composed of Fourier modes (de- 
grees of freedom). Some of these Fourier modes (those 
corresponding to the largest length scale in the system) 
appear to spend long periods of time in restricted regions 
of phase space (i.e., in some interpretations, thus break- 
ing the ergodicity assumption in statistical mechanical 
systems). As it will be clear from the present study, this 
effect is essentially the same as that which produces the 
1// spectra. 

Additionally, we make a further connection with an- 
other phenomenon observed on long times scales - the 
spontaneous flipping or reversal of some large scale fields. 
A prominent example is the reversal of magnetic field or 
its magnetic dipole in MHD systems, an effect well known 
in the context of geomagnetic studies [l6l - [l8j . 

The organization of the paper is as follows. In section 
II the model equations, the standard ideal equations of 
HD and MHD fiows, are introduced, and the numerical 
method to solve them is described. In section III the re- 
sults arc presented, divided into subsections for different 
kinds of systems. In section IV a discussion is developed, 
followed by the conclusions and summary in section V. 

II. MODEL EQUATIONS 

In this paper we consider several systems of equa- 
tions, including HD fiows with and without rotation, and 
MHD flows with and without externally imposed mag- 
netic flelds. We also consider an approximation of the 
MHD equations in the limit of strong imposed magnetic 
fields, the so-called reduced MHD (RMHD) equations. 
All systems are three-dimensional, and the fiows are in- 
compressible and ideal (zero dissipation coefficients). 

In the most general case, the incompressible HD and 
MHD equations can be written in dimcnsionless form as 

+ w X u + 2rj X u = — VT'+jxB, (1) 
at p 

— = V X (u X B), (2) 

where B is the magnetic field, u the velocity field, j = 
V X B the current density, w = V x u the vorticity, and V 
the pressure. A term considering rigid rotation with an- 
gular velocity VL can be included in the velocity equation. 
This term corresponds to the Coriolis force, with the cen- 
trifugal force absorbed into the total pressure term. The 



pressure can be obtained taking the divergence of the 
velocity equation, using the incompressibility condition 
V-u = 0, and solving the resulting Poisson equation. The 
solenoidal (V • B = 0) magnetic field includes a uniform 
part Bq and a fiuctuating part b, so that B = b + Bq. 

With Bo = 0, these equations reduce to the MHD 
equations without an external field. When B = 0, the 
equations are the ideal incompressible three-dimensional 
hydrodynamic equations, that is, the Euler equations, 
which we also consider as a case study. Finally, when 
Jl 7^ any of these systems is written in a rotating frame, 
while r2 = corresponds to the non-rotating case. 

The RMHD equations can be derived from Eqs. ([T]) 
and ([2]) for = and in the limit of strong Bq = 
-Bqz, assuming low frequencies and weak spatial gradi- 
ents along the direction of the background magnetic field 
[l9l[20l|. The equations involve potentials a(a;, y, z, t) and 
^{x, y, z, t) such that in rectangular (x, y, z) coordinates 
one has b = Vj_ x za and u ~ Vj_ x iij}, where ~ 
{dx,dy,0). To get slow dynamics as Bq ^ 1/e — > oo, for 
small e, one is forced to an ordering such that dz = 0{e). 
In this sense the RMHD equations are "quasi-two dimen- 
sional" in X and y. The dynamical equations become 
simply those of the potentials, 

1^ + u • v^uj = h ■ v^j + Boff + ^ly^uJ 

|^ + u.V^a= Bo^+^lV^a (3) 

where the electric current density is j = — V^a and the 
vorticity is a; = — 

In the subsequent section we will describe numerical 
results obtained with two different type of codes, which 
we now describe briefly. The parameters of all runs dis- 
cussed below are given in Table I for reference. 

For one set of numerical simulations, we assume pe- 
riodic boundary conditions in a cubic box, of side 2ttLo 
in each cartesian direction, with Lq the arbitrary unit of 
length. Fields can then be decomposed in Fourier modes 

u = ^Uke^'^■^ b^^bkc"'- , (4) 

k k 

where Uk, bk are the Fourier coefficients of the expansion, 
and k is the wavevector, having integer components in 
the dimcnsionless case. 

We employ a pseudospectral code to accurately nu- 
merically solve these equations. With the pseudospec- 
tral method (using full de-aliasing with the 2/3 rule) 
any quadratic invariant (like the total energy) is exactly 
maintained, except for machine round-off errors and time 
integration discretization errors. Time integration here 
is done with a second-order Runge-Kutta method, with 
a very small time step dt, to control the discretization 
error over the long simulations we carry out. Typically 
we use dt ^ 5 X 10~* which is much less than the global 
large scale turnover time Lq/u, where u is the root mean 
square (rms) velocity. As an example, in an integration 
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of 1000 unit times duration, for initial primitive fluctu- 
ation fields u and b with rms values of 1, these values 
remain equal to unity with an error less than 5 x 10~^ 
at the end of the integration. In dissipative MHD the 
energy E balance equation satisfies 



!^ = -2u{<,f> + <ca'>) (5) 

The decay in energy due to time discretization errors 
can be interpreted as numerical dissipation and for the 
example mentioned, using a time averaged value for < 

> + < ■up' >^ 10, the numerical viscosity is estimated 
as 1/ ~ 10^1". 

In a second set of numerical simulations, we consider 
spherical geometry: the HD or MHD equations are solved 
inside a sphere of unit dimcnsionlcss radius, with vanish- 
ing velocity and magnetic field at the sphere boundary. 
For this geometry, a fully spectral Galerkin code is used, 
based on a Chandrasekar-KendaUfC-K) decomposition 
of the fields. The C-K functions [Hi, [13 are 

Jj = AV X n/>j + V X (V X ril)i) , (6) 

where we work with a set of spherical orthonormal unit 
vectors {f,9,(j)), and the scalar function ipi is a solution 
of the Helmholtz equation, (\7^ + X'^)tpi = 0. The explicit 
form of tpi is 

Mr,0,<P) = C,i ji{\X,i\r)Yimi0,(t>), (7) 

where ji{\Xqi\r) is the order-Z spherical Bessel function of 
the first kind, {Xqi} are the roots of ji indexed by q (so 
that the function vanishes at r = 1), and Yim{d-,4>) is a 
spherical harmonic in the polar angle and the azimuthal 
angle (j). The sub-index i is a shorthand notation for the 
three indices {q, I, m); g = 1, 2, 3, . . . corresponds to the 
positive values of A, and q = —1, —2, —3, . . . indexes the 
negative values; finally / = 1,2,3,..., and ~l < m < I. 
The C-K functions satisfy 

V X J, = A,J, . (8) 

With the proper normalization constants, they are an 
orthonormal set that has been shown to be complete [1^ . 
The values of |Ai| play a role similar to the wavenumber k 
in the Fourier expansion. Note that boundary conditions, 
as well as the Galerkin method to solve the equations 
inside the sphere using this base, were chosen to ensure 
conservation of all quadratic invariants of the systems, 
crucial for our present study of ideal flows for long times. 
More details about the technique to numerically solve the 
HD and MHD equations in this spherical geometry can 
be found in [23,[2|. 



III. RESULTS 

A. Three-dimensional MHD in a box and in the 
sphere 

For three-dimensional incompressible ideal MHD with 
no mean magnetic field, there are three quadratic invari- 
ants: the total (kinetic plus magnetic) energy per unit 
mass 

E = ^{\u\^ + m = E., + E,, (9) 

(with (...) denoting a spatial average), the cross helicity 

= (u • b) , (10) 

and the magnetic helicity 

i7,„ = (a-b) (11) 

, where a is the vector potential such that V x a = b. The 
robustness of the Gibbs equilibrium ensemble predictions 
for this system is well established [26l - [30j . The expecta- 
tion value of spectra are readily obtained in the Gibbs 
ensemble using Lagrange multipliers associated with each 
conserved quantity. These expectations are well verified 
with numerical simulations [29l . |30| . Usually the Gibbs 
ensemble has been viewed as a predictor of the direction 
of spectral transfer, although more recently it has also 
been found that the ideal Gibbs-Galerkin system shares 
additional characteristics with the dissipative turbulent 
system at short and intermediate time scales and length 
scales. One of the characteristics found is that even in the 
ideal truncated case, cascades develop during transients 
f [3ll - [33 |l. For long times the system goes to solutions 
with zero flux but any perturbation of the system away 
from the zero flux solutions (e.g., by thermal fluctua- 
tions) is corrected by transient non-zero fluxes associated 
with the non-linear interactions ([39|). 

For wavenumber fc >> 1 (the wavenumber k in the fol- 
lowing discussion should be considered equivalent to |Ai| 
for the spherical case), the Gibbs ensemble predicts in 
three-dimensions an omnidirectional spectrum going like 
k^. Attaining this equilibrium prediction is frequently 
called "thermalization" of the large wavenumber modes, 
as this corresponds to equipartition of energy among all 
individual modes, i.e., a flat modal spectrum. At the fun- 
damental modes k = 1 (the largest possible wavelength 
modes in a finite size system), condensation is predicted, 
according to the established value of Hm- When Hm ^ 0, 
condensation at the fc = 1 mode occurs, and this has been 
the base for prediction [2^ of an inverse cascade of mag- 
netic helicity in dissipative MHD (i.e., when dissipation 
and forcing is added to the ideal MHD equations). 

Furthermore, when both Hm ^ and He ^ 0, the 
condensation of magnetic helicity induces a partial con- 
densation of the cross helicity [30] ■ For such cases the 
largest scales are expected to contain signatures both of 
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magnetic helicity (helical b) and of Alfvenic correlation 
(u oc b). 

We note here that the ideal model can be viewed as a 
dynamical model of the nonlinearities that drive turbu- 
lence. It is a simplified model because it becomes Gaus- 
sian, and also it lacks a preferred average direction of 
spectral transfer. But it is worth remembering that in 
"real turbulence" , with dissipation, there are large num- 
bers of couplings that take energy to higher k, and also 
large numbers of couplings that take energy to lower k. 
These two types of couplings are almost in balance. How- 
ever, transfer to higher k dominates slightly - this is just 
what the direct cascade is (an inverse cascade would be 
the case where the transfer to lower k dominates slightly). 
This has been observed in numerical simulations, as well 
as in laboratory experiments. In this regard, the ideal 
model is not "real turbulence" but it shares some of its 
properties ([HI). 

An assumption for the equilibrium ensemble predic- 
tions is the property of ergodicity. The required property 
is that, in a long period of time, the system, defined by 
the set of real and imaginary parts of its Fourier coeffi- 
cients (i.e., the dynamical degrees of freedom of the ideal 
MHD system) will visit all accessible regions in complex 
phase space. A point in phase space is "accessible" if it is 
permitted by the values of the invariant quantities. This 
is well verified for k > 1 modes, which evolve on relatively 
fast time scales. However, for the longest wavelength 
k ~ 1 modes there is an apparent breaking of ergodicity 
that has been identified in simulations [iJ: Eli- There 
is however evidence that ergodicity is restored at very 
long times [ij]. As a result, the fc = 1 modes seem to 
wander for very long periods of time in restricted regions 
of phase space, until at some point, a possibly sudden 
"hopping" [l3| occurs and a new period of wandering oc- 
curs in another restricted region of phase space. As has 
been recently also pointed out the time duration of 
these wandering periods is related to the dimension of the 
system (i.e., the finite number of modes assumed for the 
numerical simulation). This is to be expected since the 
delayed ergodicity is related to condensation, and for any 
system with the same values of the ideal invariants, the 
condensation is more complete for systems with larger 
numbers of degrees of freedom [30| . 

We will focus here then on the time behavior of a sin- 
gle fc = 1 mode, and the point we make is that there 
is a connection between this aperiodic cycle of wander- 
ing and hopping behavior, and the 1/f power frequency 
spectrum already observed in dissipative MHD The 
time behavior of the fc = 1 mode can be interpreted 
as the time behavior of the large scale magnetic field, 
i.e., it corresponds to the time behavior of the fiuctuat- 
ing magnetic field after a filtering of the smaller scales 
(with fc > 1) is performed. This connection between 
fc = 1 and the apparent or observable large scale mag- 
netic field becomes sharper as the condensation becomes 
more complete. In dissipative MHD, this filtering pro- 
cess occurs naturally through dissipation which damps 



the small scales. In ideal MHD, however, the thermaliza- 
tion of the small scales, with a spectrum going like fc^, 
tends to obscure the time behavior of the fc = 1 mode if 
a power frequency analysis is applied to the full fluctuat- 
ing magnetic field. As a result, we shall look at the time 
behavior of a single fc = 1 mode in this case, in order to 
properly see the long time fluctuations that develop. In 
what follows, we thus present a series of results for the 
time behavior and the frequency spectrum of the fc = 1 
mode for different numerical simulations of ideal MHD, 
with first a fixed system size (fixed N, with = to- 
tal number of modes) and varying the magnetic helicity, 
and secondly by varying the number of modes but 
with approximately fixed magnetic helicity. Our work- 
ing hypothesis is that these scalings can be understood 
in terms of scalings within the Gibbs ensemble [30| : For 
fixed Hm/E, increasing intensifies the condensation 
until all Hm asymptotically resides in fc = 1. For fixed 
N, increasing Hm/E increases condensation until as this 
ratio approaches its maximum value, all excitations con- 
dense to fc = 1. 



1. Periodic box 

We begin by considering the behavior in time of a large 
scale mode in a sequence of three incompressible ideal 3D 
MHD simulations with increasing These are fur- 

ther described in Table I, and have H„i = -0.008, 0.129, 
and —0.395, a fixed simulation size of = 16'^, and 
fixed total energy equal to 1. The fluctuations are ini- 
tially equipartitioned £"„ ~ Eh, and concentrated in a 
range of wavenumbers 1 < fc < 4. In particular Fig- 
ure H]) illustrates the temporal behavior of the fc = 1 
modes, choosing the real part of ^^(k = 1, 0, 0) as an in- 
dicator of the behavior of the large-scale modes. Indeed, 
the same behavior is obtained for other components, bx, 
by, or for other directions in k. except of course for di- 
rections for which k is parallel to the field component, 
which are identically zero due to the V • b = condition. 
The different panels in Figure ([T]) correspond to increas- 
ing values of |i?rn|- It is apparent that as the magnitude 
of magnetic helicity is increased, long period fluctuations 
(as long as 1000 unit times) begin to appear. Such fluctu- 
ations are not observed in the time series of modes with 
larger wavenumber (i.e., smaller spatial scales), and are 
characteristic of the largest scale Fourier modes in the 
system. Note that the value of the magnetic field fluctu- 
ation amplitude increases with Hm, consistent with the 
condensation phenomenon. 

Qualitatively, it would be reasonable to say that the 
low helicity case in Figure ([T|) (top) is more "station- 
ary" , as the fluctuations in that case approach a zero 
mean within say 100 time units or less. On the other 
hand the two higher helicity cases exhibit coherent fluc- 
tuations at a scale of hundreds, or even thousands, of 
time units. The suggestion is that low frequency oscilla- 
tions are becoming more dominant with higher helicity. 
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Indeed, there is clear evidence of this in the frequency 
spectra, shown in Figure ([2]) for each of the three time 
series given in Figure ([T]). The increasing power in the 
low frequency part, and the emergence of a ~ 1// power 
law at the very low frequencies << 1 , are associated 
with increasing \Hm\. We note here that often in the lit- 
erature, 1// is loosely used to refer to any spectrum of 
the form /~", with < a < 2 (i.e. omitting both white 
noise and Brownian motion). 

The important point that we want to stress here is not 
the exact power law index that fits the spectra (which 
in fact is dependent on the time duration of the sim- 
ulations, because when longer time fluctuations appear 
a more extended run would be needed), but rather the 
fact that there is no obvious reason for which fluctuations 
with time scales orders of magnitude longer than the unit 
time should appear here. The longest time scale based 
on local time arguments is = L/ul, which for the 
largest length scale oi L ~ 2tt (box size) corresponding 
to fc = 1 and for = 1, is Tl ~ 2tt. Frequencies below 
Jl = = l/(27r) should normally have a flat power 

spectrum (white noise), indicating the existence of a de- 
fined correlation time. However as the Figure ^ shows, 
for the cases where long time fluctuations appear, the 
spectrum to the left of the /l = l/(27r) is far from being 
flat. This is indicated in the panels of the figure with 
a vertical line at the frequency 1 / (27r) and a horizontal 
line at the corresponding value of the power spectrum 
for that frequency. We identify the difference between 
the observed power spectrum and a flat power spectrum 
as the range oil/ f spectra in each case (again, using a 
loose definition for 1//). 

This range of the 1// spectra is short for the lower 
values of Hm, with corresponding time periods (T ^ 1//) 
of the order of 10 unit times, and increasing to T w 
100 - 500 for H„, = 0.129, and as long as T « 500 - 1000 
for = 0.395. 

It is interesting to observe the autocorrelation func- 
tion for each of the time series in the cases shown, as a 
complementary way to see the appearance of long time 
fluctuations. Here the autocorrelation function is defined 
as 

C(<)=<fe(to)5(io+t)> (12) 

where b represents a cartesian component of a magnetic 
field mode (for instance the z component of the real part 
of the magnetic field Fourier mode for k = (1, 0, 0)), with 
subtracted mean and normalized so that < 6^ >= 1, to is 
an arbitrary origin in time, t is the time lag and < ... > 
denotes a time average. 

This is shown in the panels of Figure ([3]). The case 
with smallest presents a compact correlation func- 
tion, localized within about 10 unit times, whereas the 
cases with larger H„i show a much broader correlation 
function. A correlation time can be obtained as 

tc= / C{t)dt (13) 
Jo 



where Tf is the final time of the run, assuming it is a 
large time. 

This quantity is well defined if the correlation function 
is confined within a finite range of time lags, decaying 
faster for long times. However, if the correlation func- 
tion does not decay faster, and in depends on the time 
duration of the series, then the correlation time is not well 
defined. In particular, for an exact 1// power spectrum 
there is no single correlation time that can be defined. 
The values indicated in the plot for the three Hm cases 
of tc = 14, 189, 346 show again that longer time fluctua- 
tions (and correlations) appear as Hm is increased. 

Recalling that condensation in the ideal Gibbs ensem- 
ble intensifies as the number of degrees of freedom (~ N^) 
increases, we now explore whether the emergence of long 
time scale coherent fiuctuations and the associated 1// 
spectra behave the same way. The next series of plots 
in Figure ^ show the time behavior of a magnetic field 
component of a single mode k = (1,0,0) varying the size 
of the simulation of 3D MHD, for = 16^,323, and 
64^, and for low values of Hm = 0.027, 0.027, and 0.015, 
respectively, for each N . These results show that long 
time fluctuations appear even for low values of Hm when 
the size of the system is increased. We note that for the 
= 64 case a longer time range is needed to observe 
fluctuations to average to zero. This is indeed observed 
in this run when it is extended until t > 10000 but not 
shown here for consistency with the shorter time range 
selected for the rest of the cases in this Figure. Figure ([S]) 
shows the corresponding power frequency spectra P{f) 
for each size N. The range of low frequency 1// noise in- 
creases from T - 10 for AT = 16 to r - 2000 for A^ = 64 
(compare Figures ^ and ([5])). 

Furthermore. Figure ([6]) shows the autocorrelation 
function and the obtained correlation time tc for each 
time series for different values of A^. Also evident in this 
figure is the appearance of long time fluctuations as A^ is 
increased through the broadening of the autocorrelation 
function. 

Another view of the long time fluctuations can be 
obtained from the time evolution of individual Fourier 
modes in a phase space plot in the complex plane [14i.il5|. 
For instance, the 6^(1,0,0) mode for the 32'^ run in the 
Hm — 0.027 case is shown in Figure ([7]). For comparison, 
the behavior in the complex plane of the 6z(2, 0, 0) mode, 
with a larger fc, for the same run is shown in Figure ([5|). 
The long time fluctuations of the 6z(l,0, 0) mode corre- 
spond to long periods of time spent in a restricted region 
in the complex plane, as contrasted by the quick filling of 
the complex plane allowed region for the 6z(2, 0, 0) case. 
This phenomenon has been called [l^ "delayed ergodic- 
ity," because the ergodicity property of the fc = 1 seems 
to be broken only temporarily, as the mode spends long 
times in a region of the complex plane thus not filling 
the entire space. For longer times, the "hopping" of the 
mode between different regions starts filling the space. 
As observed here, this corresponds to long time fluctua- 
tions of the time series, and to the appearance of the 1// 
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FIG. 1: Time series of a component of the magnetic field for 
the k = (1,0,0) mode in a MHD run with 16"^ modes and 
with different values of the magnetic helicity: H,n = —0.008 
(top), 0.129 (middle), and -0.395 (bottom). 



power law in the power frequency spectrum. 

As can be seen in the time series plots, the compo- 
nent of the large scale magnetic field shown, which is a 
dominant contribution to the global magnetic field when 
condensation is strong, progresses for long periods of time 
without changing sign, and then experiences a reversal in 
sign, followed by another long period of time without sign 
change (see e.g.. Figures ([T]) and (|4])). We will come back 
to this "reversal" phenomenon in spherical runs, but it 
can already be seen that this is part of the same long 
time fluctuations phenomenon reported as 1// noise. 

Another interesting diagnostic is shown in Figure 
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FIG. 2: Power frequency spectra of the time series in Fig- 
ure with different values of the magnetic helicity Hm. = 
—0.008,0.129, and —0.395. The reference vertical line corre- 
sponds to a frequency value /l ~ l/(27r) (see text) and the 
horizontal line to the value of P^Jl)- 



where the real and imaginary parts of the complex ampli- 
tudes of several field components of the mode k = (1, 0, 0) 
are shown as a function of time, for the same run (A'^^ = 
32'^ and Hm = 0.027). This is a case with long time 
fluctuations and delayed ergodicity. Apparently, from 
our analysis the fc = 1 magnetic modes satisfy particu- 
lar equilibrium configurations, namely the field is quasi 
force-free. A force-free magnetic field satisfies j x b = 0, 
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FIG. 3: Self-correlation function for the time series of mag- 
netic field in Figure ([1]), with different values of the magnetic 
hehcity = -0.008, 0.129, and -0.395. A correlation time 
tc (see text) is indicated in each case. 



FIG. 4: Time series of a component of the magnetic field 
for k = (1,0,0) mode in the MHD runs with different reso- 
lution iV^ = 16^,32^, and 64^ and magnetic helicity Hm = 
0.027, 0.027, and 0.015. 



that is, the Lorcntz force term in the momentum equation 
is zero. If the velocity field u is zero, this means that for 
an ideal flow (no viscosity or diffusivity) the magnetic 
field will remain force-free in time, and all non-linear 
terms will be zero. Force-free states also correspond to 
maximum allowable values of the magnetic helicity (see 

(H). 

In terms of Fourier modes, if a single k mode is in a 
force-free state, i.e., jk x bk = 0, then, since jk = ik x bk. 



it must satisfy ik x bk = Abk which implies that ~ k'^, 
that is, A — ±fc. For the fc = 1 modes this means 
A = ±1 and, since these modes are equal to the Carte- 
sian versors (unit vectors), it also implies some relations 
between the components of bk- Taking, for example, 
k = (1,0,0), these relations are Im{bz) = ±Re{by) and 
Re.{bz) = T^w(6y), where the field components are in 
k-space (notice that for this mode, hx — Q because of 
the divergence-free condition.) This is immediately seen 
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FIG. 5: Power frequency spectra of the time series in Figure 
(Q , with different resolution iV^ = f6^32^ and 64^ and mae- 
netic helicity H™. = 0.027, 0.027, and O.Of 5. The reference 
vertical line corresponds to a frequency value /l = l/(27r) 
(see text) and the horizontal line to the value of 
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FIG. 6: Self-correlation function for the time series of mag- 
netic field in Figure Q, with different resolution A*"^ — 
16^,32^, and 64^ and magnetic helicity = 0.027,0.027, 
and 0.015. A correlation time tc (see text) is indicated in 
each case. 



to be the condition for a circularly polarized fluctuation 
in a k-aligned coordinate system - that is, if we let 
k = (1,0,0) in a right-handed Cartesian system, then 
the above conditions are equivalent to bk = e'^'=(0, 1, ±z) 
for arbitrary phase (f>. The plots in Figure ([9]) shovif that 
the imaginary and real parts of these field components 
exhibit this correlation - they are very highly correlated 
in time (light gray and black are used for each corre- 



sponding component and arc almost indistinguishable). 
In fact, values of the correlation coefficient > 0.95 are 
obtained in each case, corresponding to near maximum 
helicity. 

Similar relations can be found for the k = (0, 1, 0) and 
k = (0, 0, 1) modes, and strong correlations > 0.95 are 
also found for the imaginary and real part of the corre- 
sponding field components in time (not shown). This is 
true for all the previous runs for which long time flue- 



9 



t=40-125 



t=40-250 



-0.05 0.00 0.05 

Rc|b,[k=(l,0,0)]| 



(=40-375 




-0.05 0.00 0.05 

Rc|b,[k=(l,0,0)]| 




-0.05 0.00 0.05 

Re|b,[k=(l,0,0)]l 



t=40-500 




-0.05 0.00 0.05 

Re|b,[k=(l,0,0)]l 



FIG. 7: Complex phase space trajectory for the k = (1, 0, 0) 
mode in the MHD 32^ rmi, with Hm = 0.027, for different 
intervals of time. 
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FIG. 8: Complex phase space trajectory for the k = (2, 0, 0) 
mode in the MHD 32'^ run, with magnetic helicity Hm ~ 
0.027. Notice the noticeably larger wandering in phase space 
when compared to Figure ((TJ, indicative of clearer ergodicity. 



FIG. 9: Top: real part of the by component of the k = (1, 0, 0) 
mode vs time (black), and imaginary part of bz for the k = 
(1,0,0) mode vs time (light gray) for the MHD 32^ run with 
magnetic helicity Hm = 0.027. Bottom: imaginary part of 
—by for the k = (1,0,0) mode vs time (black), and real part 
of bz for the k — (1, 0, 0) mode vs time (light gray). 



tuations of the k = 1 modes are observed. However the 
circular polarization is not found for larger k modes (e.g., 
for k (2, 0, 0)). This is a very special property of the 
k = 1 modes, that they evolve in a quasi force- free state 
in time. This condition is consistent with the ensemble 
average predictions of the Gibbsian theory [2^, [sO] , but is 
not imposed by it as an exact condition due to allowance 
for fluctuations about the equilibrium expectation. We 
will come back to this issue in the discussion section. 



2. Spherical geometry 

The next example corresponds to a series of spherical 
MHD runs. Figure PH)) shows three runs with a fixed 
value of H„i = 0.03, and with three different values of 
Qmax = 5, 7, and 8, corresponding to increasing number 
of degrees of freedom in the model. Here, Qmax is the 
maximum value of the q index corresponding to the radial 
number in a spherical harmonic expansion [see Eqs. 1^- 
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([7]) and text therein] ; the maximum values of I and m are 
also increased accordingly. In this particular geometry, 
the maximum possible value of helicity for a flow with 
unit energy is ~ 0.22. Again, as in the periodic box runs, 
it is seen that increasing the size of the system shows the 
appearance of long time fluctuations and a corresponding 
range of 1// power frequency spectrum, as shown in the 
plots of Figure ([TT|) . 

Here we show again in Figure (jll[) references vertical 
line for the frequency = 1/2 based on the largest local 
non-linear time T/, = L/Ul = 2 which corresponds to 
structures with the diameter of the sphere L = 2 (unit 
radius) and a unit root mean square velocity [/^ = 1. 

Another effect is shown in Figure ([12]) which corre- 
sponds to the time behavior of a single mode, by in- 
cluding non-zero rotation, with f2 — 16. Results for size 
Qmax = 5 with $7 = and O = 16 arc shown. It is seen 
that longer time fluctuations appear with the addition 
of rotation in the system. The corresponding frequency 
spectra in Figure (|13p also show a wider range of 1/f 
noise for the rotating case. 

The long time fluctuations for the case of qmax — 5 and 
Hm, = 0.03 with rotation (sec Figure ([T^ ) appear as sign 
reversals of the z-component of the magnetic field, with 
excursions with periods of the order T ^ 50 — 200. In 
fact, another quantity that can be studied for this system 
is the magnetic dipole, which is essentially a weighted 
moment average of the magnetic field, which tends to 
highlight the low q modes. We defer a detailed analysis 
of the magnetic dipole reversals to another paper, but 
we point out that this long time fluctuation phenomenon 
is directly related to 1/f noise and delayed ergodicity of 
the q ~ I mode in complex phase space. As an example, 
a phase space plot of the temporal behavior of the q = 1 
mode is shown in Figure It can be seen that this 

mode spends a long time in a restricted region of phase 
space, thus delaying overall ergodicity. 



B. Three-dimensional MHD and RMHD with a 
background magnetic field in a cubic box 
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1. MHD with a background magnetic field 

In the presence of a background uniform magnetic field 
Bo, the 3D MHD equations lose one of the quadratic 
invariants, the magnetic helicity, so the two invariants 
that remain are the total energy (magnetic plus kinetic) 
and the cross helicity. As will be discussed in a following 
section, this fact can have an influence on the interaction 
of the lowest k = 1 mode with the remaining modes in the 
system, and correspondingly on the 1/f power spectrum. 

Results for simulation runs of ideal 3D MHD with a 
background magnetic field in the y-direction Bq^ = 8 

1/2 

(as compared to fiuctuations r.m.s. values of (6^) = 
brms = 1) are shown in Figure ([TSl) . for the time evolution 
of a fluctuating component of the magnetic field of the 
k — (1,0,0) mode, and for the corresponding frequency 



FIG. 10: Time series of the bz component of the magnetic 
field for the {q,l,m) = (1,1,0) mode in the spherical MHD 
run with qmax = 5, 7, and 8, and magnetic helicity Hm = 0.03. 



spectra in Figure (jl6p . This corresponds to a simulation 
of size 32"^. Long time fluctuations of the order of 1000 
unit times appear, and a corresponding strong enhance- 
ment of very low frequency power, associated with a 1/f 
spectrum, is clearly observed. This case of MHD with a 
background magnetic field has been thoroughly studied 
for the driven/dissipative case [l[. Here, it can be seen 
that the 1/f spectrum appears even in the ideal case 
and we conclude that it is a property of the non-linear 
interactions among the modes in the system, and docs 
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FIG. 11: Power frequency spectra of the time series in Figure 
(fTO)l . for q,nax = 5,7, and 8, and Hm = 0.03. The reference 
vertical hne corresponds to a frequency value /l = 1/2 (see 
text) and the horizontal line to the value of P(/_l). 



not depend on the presence of dissipation in an essential 
way. The complex plane phase space trajectory for this 
case is shown in Figure ([TT]) . and delayed ergodicity of 
the k = (1, 0, 0) mode is also clearly observed. 




FIG. 12: Time series of the bz component of the magnetic 
field for the {q,l,m) = (1,1,0) mode in the spherical MHD 
run with zero rotation (fi = 0, top), and 57 = 16 (bottom) 

and Qrnax ~ 5. 



2. Reduced MHD 

Another related case of interest is given by the RMHD 
system, which is, as discussed in Section II, an approxi- 
mation to the low-frequency dynamics of the MHD equa- 
tions with a large background magnetic field. Results for 
a RMHD case are shown in Figure (|18p . namely time 
series for a component of the fluctuating magnetic field 
for the k = (0, 1, 0) mode, and corresponding power fre- 
quency spectrum in Figure (jl9p . This is for a case with 
a background magnetic field in the a;-direction Bq^ = 8 
and size 32'^. As can be seen, long time fluctuations are 
evident, with periods on the order of 1000 unit times. 
The RMHD system admits two known quadratic ideal 
invariants, the total energy and the cross hclicity. It 
is interesting however that an additional quantity, the 
mean square vector potential, which is strictly invariant 
in ideal 2D MHD, behaves as a quasi- invariant in RMHD 
[3^. By this we mean that this quantity remains statis- 
tically constant for long periods of time. The connection 
of this with the emergence of long time fluctuations will 
be further discussed in section IV. 
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FIG. 13: Power frequency spectra of the time series in Figure 
p2|l . with rotation f2 = and 16 and qmax = 5. The reference 
vertical hue corresponds to a frequency value /l — 1/2 (see 
text) and the horizontal line to the value of P{fL)- 



C. Three-dimensional hydrodynamic with and 
without rotation 

There are two quadratic invariants in ideal three- 
dimensional hydrodynamics, the kinetic energy 



Eu^l{\u\'), 



and the kinetic helicity 



(14) 



(15) 



These are also invariants of the 3D HD equations with 
addition of constant rotation. 

In the absence of rotation, these invariants do not con- 
dense to the longest wavelength, or engage preferentially 
in back-transfer to the lowest k values [32,, jJTj. Lack- 
ing these tendencies a rationale is lacking for expect- 
ing an inverse cascade to large scales in the dissipative 
case. In the presence of rotation, however, and in the 
ideal case, there is a transient reduction to quasi two- 
dimensional behavior and a transient condensate, asso- 
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FIG. 14: Complex phase space trajectory for the (g, I, m) = 
(1,1,0) mode in the spherical MHD run, with Qmax = 5 and 
rotation Q — 16. 
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FIG. 15: Time series of the z-component of the magnetic 
field for the k = (1,0,0) mode in the MHD 32^ run with 
background magnetic field Bq^ = 8 and brms = 1- 



ciated with the quasi-conservation of the energy in two- 
dimensional modes (modes with — 0, where parallel 
refers to the direction of the rotation axis) [H, HI] . These 
modes correspond to the so-called slow manifold of the 
system. In the forced-dissipative rotating case, resonant 
interactions transfer energy toward these modes, also re- 
sulting in two-dimensionalization and in the development 
of an inverse energy cascade [13, 113 ■ 

The ideal behavior discussed above is similar to the 
behavior of MHD in the presence of a background mag- 
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FIG. 16: Power frequency spectra of the time series shown 
in Figure (|15|l . The reference vertical line corresponds to a 
frequency value /l = 1/ {2tt) and the horizontal line to the 
value of P(/l). 



FIG. 18: Time series of the z-component of the magnetic 
field for the k = (0, 1,0) mode in the RMHD 32^ run with 
background magnetic field Bq^ = 8. 
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FIG. 19: Power frequency spectra of the time series shown 
in Figure p8p . The reference vertical line corresponds to a 
frequency value /l — l/(27r) and the horizontal line to the 
value of P(/i). 



FIG. 17: Complex phase space trajectory for the k = (1, 0, 0) 
mode in the MHD 32'' run with background magnetic field 
Bo„ = 8. 



netic field (see previous subsection), which also reduces to 
a quasi two-dimensional dynamical behavior and has the 
mean square vector potential as a quasi- invariant [28l[36l |. 
As will be discussed further in the next section, this has 
an effect on the emergence of long time fluctuations (see 
also Hi). 

The time evolution of a velocity component for the k = 
(1, 0, 0) mode with no rotation and with rotation fi, = 16 
(along the z-axis) in a periodic box is shown in Figure 
(j20p . and the corresponding frequency spectra are shown 
in Figure (|2ip . Long time fluctuations are much more 
apparent in the case with rotation, where fluctuations 



arc observed with periods of the order of T 50 — 100, 
together with a corresponding enhancements in the low 
frequency l/f range of the spectrum. Similar results 
are obtained for simulations in spherical geometry (not 
shown) . 



IV. DISCUSSION 

All the simulations discussed in the previous section 
are summarized in Table I. We summarize now the results 
and discuss the related phenomena of emergence of long 
time fluctuations, delayed ergodicity, and corresponding 
l/f power law in the frequency spectra. 

The flrst examined case corresponds to ideal 3D MHD. 
This system shows long term memory, l/f noise, and de- 
layed ergodicity in the k = 1 modes. As pointed out. the 
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TABLE I: Summary of all runs, with A'^'' the total number of grid points, Hm, and Hv the magnetic and kinetic helicity, and 
n or Bo the respective amplitude of the imposed solid body rotation or uniform magnetic field. The first column indicates 
what equations have been integrated, SMHD being the MHD equations in spherical geometry with qmax corresponding to the 
truncation of the C-K functions (see equation (5) and text); all other runs are computed in periodic cubic geometry. "ROT" 
stands for runs with rotation, and Bq stands for runs with an imposed magnetic field. The last column gives the estimated 
strength and range of the observed 1// noise spectrum (see Figures). 
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FIG. 20: Time series of the ^-component of the velocity field 
for the k = (1,0,0) mode in the HD 32'' run with rotation 
n = (top) and Vly = 16 (bottom). 



FIG. 21: Power frequency spectra of the time series of Figure 
(|20p . The reference vertical line corresponds to a frequency 
value /l = 1/ (27r) and the horizontal line to the value of 
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system has three quadratic invariants, the total energy, 
cross hehcity, and magnetic helicity, and in a statistical 
steady state the amplitudes of Fourier modes arc con- 
trolled by the Gibbs ensemble prediction. Specifically, 
the magnetic helicity allows condensation at the lowest 
wavenumber mode. This happens for non-zero values of 
and it becomes more intense as the absolute value 
of Hrn is increased. Condensation also become more in- 
tense as the number of modes ~ is increased for fixed 
Hjn/ E. The k = 1 mode is special then in the dynamics 
of this system. As shown, this mode is in a quasi force- 
free state, so its evolution is slow, weakly coupled with a 
sea of lower amplitude modes at much larger wavenum- 
bers. The force- free property is prescribed by the Gibbs 
ensemble solution; the modes with fc = 1 have maxi- 
mum helicity, the magnetic field is parallel to the current 
density, and these largest scale fluctuations are circularly 
polarized. 

The coupling between the fc = 1 modes and the small 
scale modes is defined by triads of wavenumbers, con- 
structed with a fc = 1 mode and two large wavenum- 
ber modes. The time evolution of these interactions 
is controlled by the (comparatively small) amplitude of 
the larger wavenumber modes, but the large length-scale 
(small fc) of the lowest wavenumber mode. As pointed 
out in [l|, this can be seen from the expression for this 
type of interaction which is of the (schematic) form 

^^-^kY: <i)Kp) , (16) 

k=p+q 

where b{k) , u{q) , b{p) are generic Fourier mode ampli- 
tudes, with the constraint that k = p -I- q. In partic- 
ular, we consider the lowest wavenumber mode fc = 1. 
If the triadic interaction is local, then k p ^ q hy 
definition and the timescale of that interaction is given 
by [fcu(fc = 1)]^^ ~ 1, whereas if the interaction is non- 
local, then p, (j'3>fc = l, and the timescale is 
[ku{q)b{q)/b{k = 1)]^^ which is much longer than the 
local timescale since u{q), b{q) <^ u{k = 1), b{k = 1). 

The fact that 1// noise is also observed in cases with 
low magnetic helicity but increasing size N suggests how- 
ever that even in cases with amplitude equipartition 
among modes (i.e., when magnetic helicity is small, as 
indicated by the Gibbs ensemble predictions) , the sea of 
large wavenumber modes is again slowly modifying the 
dynamics of the fc = 1 mode. This idea is supported 
by previous studies [11] in which the effect of the large 
wavenumber modes on the lower wavenumber modes is 
modeled through an effective viscosity, even though the 
systems are strictly ideal, like the Euler equations. This 
is also similar to ideas suggested by low dimensional dy- 
namical systems in connection with the reversals of the 
geomagnetic field, see e.g. [13, 13, with bi-stable states 
driven by noise (in this case, the large wavenumber modes 
would act as the driving noise for large-scale behavior). 

These results, besides being obtained for different res- 
olutions and values of the invariants in ideal cases (i.e.. 



without viscous or external forcing, indicating the long 
term behavior is intrinsic to the system of equations), 
are also obtained for two different geometries: in peri- 
odic boxes and in spheres. As a result, we conclude that 
the boundary conditions do not seem to affect the long 
term behavior. 

As mentioned before, the 3D MHD system has an ideal 
invariant that experiences a condensation to the longest 
wavelength modes of the system, and in the dissipa- 
tivc driven case is expected to be involved in an inverse 
cascade. Interestingly, another set of results observed 
here suggest that certain systems like MHD with a back- 
ground magnetic field or HD with rotation also allow for 
the emergence of long time fluctuations and 1// behav- 
ior. These systems do not have an ideal invariant that 
condenses to the lowest wavenumber mode in the Gibbs 
ensemble. There is nevertheless a slow manifold of modes 
(the two-dimensional modes) that is distinguished from 
all other degrees of freedom. We argue that this slow 
manifold is controlling the emergence of long time fluc- 
tuations in the fc = 1 modes, as observed in the results. 
It is interesting to note that these systems are charac- 
terized by the existence of a quasi- invariant. By this we 
mean a slowly varying but not strictly conserved quan- 
tity, such as the square vector potential (a^) in MHD 
with a background field, or the energy in modes with 
fc|l = in HD with rotation. The quasi-invariants are 
approximately conserved in the slow manifold, and thus 
introduce longer timescales in the dynamics by permit- 
ting transient condensation at the lowest wavenumber 
mode. 

Finally, a third category is that of systems with fiat 
frequency spectrum for small frequencies, which are as- 
sociated with negligible long time correlations. An ex- 
ample is 3D NS without rotation, or 3D MHD with no 
magnetic helicity and zero mean magnetic field. These 
systems have no invariant or quasi-invariant that conden- 
sates at large scales in the ideal case. 



V. CONCLUSION 

The emergence of long time fiuctuations and 1 // noise 
in the frequency spectrum of field variables is observed 
in systems with a quadratic invariant allowing condensa- 
tion at the lowest wavenumber mode, like the magnetic 
helicity in 3D MHD. This happens when the invariant is 
large but also when it is small, provided the number of 
modes (iV"^) is large enough. This happens indistinctly 
in geometries like a periodic box or a sphere, thus we 
argue that this is an intrinsic property of the non-linear 
couplings in the system and is not dependent on the ge- 
ometry, or other external properties of the system like 
driving, or dissipation, which are absent in the ideal sys- 
tems analyzed here. 

We can conjecture therefore that long time fluctuations 
will be also observed in ideal 2D MHD and 2D HD, where 
a quadratic invariant allowing condensation at the lowest 
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wavenumber mode also exists. A previous indication of 
this is given in [l| where the driven-dissipative case for 
these systems is studied and showed to have 1// noise. 
Other non-Unear systems that do not belong to fluid dy- 
namics, but have Gibbsian statistical condensates and 
dissipative-driven inverse cascades, may also show this 
behavior (see for instance a case in quantum optics (4ll|). 

The 1// power spectrum is also observed in systems 
with a slow manifold and quasi-invariants, like MHD with 
a background magnetic field and hydrodynamics with ro- 
tation. We also argue that other systems, like flows in 
the geostrophic approximation, with a slow manifold dy- 
namics, may show long time fluctuations as well. 

The observed occurrence of a range of 1// frequency 
spectra in interplanetary magnetic field and density fluc- 
tuations is another case of relevance to the present dis- 
cussions This signal is observed in the frequency 
range ~ 10~^— ~ 10~* Hz in the solar wind at lAU 
and beyond, but is also observed in the coronal and in 
the solar photospheric magnetic field. It is therefore pos- 
sible that the origin of this signal is either in the solar 
dynamo or in coronal dynamics, or both. Interestingly 
one might well expect slow manifold behavior in either 
of these cases, due to rotation or regional effects of mag- 
netic helicity in the dynamo, or quasi-invariance of the 
mean square potential in the corona, which is dominated 
by a strong large scale magnetic field. It is tempting to 
also associate this phenomenon to long-time memory ef- 
fects in geophysical flows, as e.g., the ocean circulation; 
indeed, the origin of multi-decadal time-scales in the cli- 
mate evolution still remains mysterious since it is not 
directly associated with a known instability and yet it is 
well observed; and reproduced, with approximate accu- 
racy, in numerical models (see e.g. p3l - l46l |). 



Finally, three-dimensional hydrodynamic flows, which 
do not have condensed invariants or quasi-invariants, do 
not develop 1/f noise in the ideal case (see also [ij for 
similar results in the forced-dissipative case). The devel- 
opment of 1/f noise in such systems, if it happens, may 
be associated with forcing or boundary conditions. 

The main result that can be concluded here is that in 
many cases, long-time fiuctuations are intrinsically given 
by the non-linear dynamics of the system, and not con- 
trolled by external properties or dissipation. Further- 
more these studies serve to strengthen the growing un- 
derstanding that there are deep connections between con- 
densation (or quasi-invariants), and 1/f signals at low 
frequency, as well between the associated slow manifolds 
and irregular large scale stochastic reversals (or delayed 
ergodicity). A number of phenomena observed in many 
types of flows can be all understood in terms of these 
couplings, and may therefore have universal properties 
linked to the presence of invariants. 
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